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METHOD FOR EQUALISING AND DEMODULATING A DATA SIGNAL 
WHICH IS TRANSMITTED VIA A TIME -VARIANT CHANNEL 



5 The invention relates to a method for equalising and 

demodulating a data signal transmitted via a time-variant 
channel to a receiver. 

Modern data transmission procedures via time-variant 
10 channels (fading channels) are susceptible to inter- 
symbol- interference ( IS I ) or inter-channel- interference 
(ICI) . Channel estimation and equalisation are therefore 
required . 

15 Conventional methods for channel estimation and 

equalisation are based upon an estimation of the channel 
impulse response as a function in the time and/or in the 
spectral domain. The channel impulse response is 
generally estimated directly using training sequences. 

20 The channel model upon which the estimation is based, can 
model either exclusively a single time function, or may 
include various paths with a different delay using the 
conventional tapped-delay model. These models, and 
therefore also the associated estimation methods, share 

25 the disadvantage that they do not take into consideration 
the geometry of the scatterers causing the distortion. 

In the context of multi-carrier methods, e.g. OFDM, 
different Doppler shifts in the individual channel paths 

30 lead to ICI, i.e. a given carrier is influenced by 

adjacent carriers. If the real channel comprises several 
paths with a different Doppler shift, a conventional 
method with direct estimation of the channel via its 
impulse response cannot determine these different Doppler 

35 shifts. Accordingly, the ICI persists, and the receiver 
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cannot achieve optimum reception and processing of the 
signal . 

Conventional understanding of the time variation of the 
5 channel is based upon the assumption, that the impulse 
response of the channel between the training sequences 
changes only slightly or in a deterministic manner, and 
that the channel estimation and tracking algorithms used 
converge adequately . 

10 

With multi-carrier methods, e.g. OFDM, it is implicitly 
assumed that the channel is constant over a OFDM block. 
For example, a method for the equalising DVB-T based on 
the assumption of constancy is described in Burow-R; 
15 Fazel-K; Hoeher-P; Klank-O; Kussmann-H; Pogrzeba-P; 

Robertson-P; Ruf-M-J "On the performance of the DVB-T 
system in mobile environments" IEEE GLOBECOM 1998. 

With very rapidly changing channels, the methods 
20 described above require a rapid sequence of training 
sequences and/or lead to a poorer convergence of the 
channel estimation. With multi-carrier methods, constancy 
over a block, as mentioned above, is no longer 
guaranteed, and the performance of the methods decline 
25 considerably. 

The object of the present invention is therefore to 
provide a method for equalising and demodulating a data 
signal transmitted via a time-variant transmission 
30 channel of this kind, which avoids the above 

disadvantages and limitations regarding the properties of 
the channel. 
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This object is achieved on the basis of a method 
according to the preamble of the independent claim by its 
characterising features. Advantageous further embodiments 
are defined in the dependent claims. 

5 

With the method according to the invention, the channel 
impulse response is no longer used for channel 
estimation. Instead, the so-called scatterer 
coefficients, that is to say, the complex-valued 

10 attenuation, the delay and the Doppler shift in the 
channel, are used. The reflections of a signal 
transmitted between a transmitter and a receiver caused 
by so-called scatterers have a causative influence on the 
quality of the transmission channel, as described for 

15 example, in the book by Raymond Steele, "Mobile Radio 
Communications", Pentech Press, London, 1992, Section 
2.3.1. Scatterers of this kind, such as buildings or 
vehicles, distort the data signal transmitted between the 
transmitter and the receiver. Scatterer coefficients in 

20 the data signal, which are attributable to the scatterer, 
can be determined in the receiver, and the distorted data 
signal can then be equalised and finally demodulated. 
According to the invention, the channel properties are 
therefore defined by these scatterer coefficients, which 

25 can be determined in a simple manner from the distorted 
data signals received on the basis of the following 
description . 

The invention will be described below in greater detail 
30 with reference to schematic drawings of exemplary 
embodiments. The drawings are as follows: 



Figure 1 
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shows the two-dimensional arrangement of the 
scatterer with the discretised Doppler 
frequencies and delays; 



5 Figure 2 shows a search tree; 

Figure 3 shows a tree derived from the search tree of 
Figure 2 taking the coding into 
consideration . 

10 

On the basis of a two-dimensional field, Figure 1 shows 
the discretisation of the Doppler frequency f d and the 
delay t in the transmission channel for various 
scatterers. This graphic representation can be directly 

15 converted into a scatterer matrix S with the scatterer 
coefficients S (m, k) , as used in the following equations 
(1) to (4). The coefficients of the matrix S represent 
the complex-valued attenuation values (amplitude and 
phase) . The quantisation in the delay direction t and in 

20 the Doppler shift direction f d depends on the channel and 
the data transmission scheme. The maximum values K for 
the discrete, standardised Doppler shift and M for the 
discrete, standardised delay, result from the physical 
parameters of the channel. As can be seen, it is 

25 advantageous and, without restriction to generality, 

useful for the quantisation in the delay direction and 
the Doppler shift direction to be equidistant in each 
case. If no physical scatterer occurs for a given entry, 
then the corresponding scatterer in the matrix is simply 

30 set to zero. 



Fi ffure^ l shows five scatterers, of which the indices 
correspond to the position in the scatterer matrix; in 
this context, the numbering begins with 1. 
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The symmetry with reference to the Doppler shift 
(positive and negative values) is not necessary a priori; 
it is dependent upon the channel. 

5 

As a result, this physical model takes into consideration 
the geometry of the channel propagation model instead of 
the pulse responses. This geometry, and therefore also 
the delay x and the Doppler shift f d allocated to the 
10 relevant scatterer, remains practically constant for 

sufficiently long periods, because the transmitter and/or 
receiver cannot move at an arbitrary velocity and/or 
cannot perform changes of movement at an arbitrary 
velocity . 

15 

By contrast, the impulse response of the channel can, in 
principle, change arbitrarily within the permitted 
physical boundaries. The discrete impulse response can be 
calculated from the complex scatterer coefficients S (m, k) 
20 to give : 

h(m,i)=-= J^mk)e J " N 

m=0 



25 In this context, K is the maximum Doppler frequency 

occurring, m is the running index for the delay and i is 
the discrete running variable for time. h(i) is the 
resulting discrete impulse response of the channel in the 
time domain. It is observed over the length N. 

30 

The time-variant continuous impulse response of the 
channel h(x,t) is physically bounded in t and f d . 
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Accordingly, for the scatterer function, S(x,f d ) as the 
Fourier transformation of h(T,t) over t and be set to 
S(T,f d )=0 for T>T max , | f d |>f d , ma x- By analogy with the sampling 
theorem, the impulse response h(x,t) can therefore be 
5 presented completely through sampled values within the 
frequency domain, so that (1) is obtained as a discrete 
presentation of the channel. 



The maximum likelihood approach for determining the 
10 scatterer-coef f icient matrix S in the time domain is 

obtained by minimising the following expression according 
to the scatterer coefficients: 



15 



N-l 



i)-*f i d(i-m)f j S(m,k) l 

\N m =0 k=-K 



J2x- 



(2) 



In this context, it is implicitly assumed, that the 
transmitted data symbols d(i-m) are known. r(i) is a 
sample of the signal received. 



20 The variables r(i) and d(i-m) are defined within the time 
domain . 



The data symbols are either assumed to be known directly 
as a training sequence or they are determined from the 
25 signal received using the method described below. 

The scatterer coefficients are preferably estimated in 
the time domain with data transmission schemes, which 
operate within the time domain. Such methods include, for 
30 example, single carrier methods with PSK or QAM 
modulation . 
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In the case of multi-carrier signals with known 
transmitted symbols, the estimation could also be carried 
out within the time domain, because the transmission 
signal is previously known. 

5 

The modulation scheme can be taken into consideration in 
equation (2), in that the data symbols d(i-m) carry the 
relevant signal form of the modulation type used, 
optionally with partial response pulse shaping. Channels 
10 with a large memory, i.e. with a long pulse duration, can 
be equalised by a corresponding choice of the maximum 
delay M. In this context, the duration of observation N 
is naturally also of a corresponding length. 

15 An estimation can be implemented in the frequency domain 
in a similar manner to equation (2). In this context, the 
following equation is obtained: 

N-\ II i K M-l _ j2mn — 1| 2 

£ \R{n) - £ £ - k)S(m, k)e (3) 

20 

The variables R(n) and D(n-k) shown in (3) are defined 
within the frequency domain. 

The scatterer coefficients are preferably estimated in 
25 the frequency domain with data transmission schemes, 

which operate within the frequency domain. Such methods 
include, for example, multi-carrier schemes such as OFDM 
with the DVB-T method. 

30 As for an estimation within the time domain, the data 

symbol D(n-k) can carry the signal form of the modulation 
type used, presented in this context, within the 
frequency domain. 
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As can be seen from equations (2) and (3), for the 
estimation of scatterer coefficients, the transmitted 
data are assumed to be known. The estimation is carried 
5 out over N samples in the time domain and/or N spectral 
components in the frequency domain. 



Normally, at the beginning of a data transmission, a 
known symbol sequence is transmitted, which is used for 

10 synchronisation. Following this, in the case of unknown 
data sequences, the receiver must track the estimation of 
the channel and/or, with a new transmission of 
synchronisation information or training symbols, re- 
estimate and/or adapt the convergence behaviour of the 

15 estimation and tracking algorithm. 



Estimation of the scatterer coefficients is preferably 
carried out by means of a recursive Kalman algorithm or 
an RLS algorithm, in which, after initialisation with the 
20 known symbol sequence, the channel is tracked with 
initially unknown sequences . An RLS algorithm for 
determining the scatterer coefficients reads, for 
example, as follows: 



25 K(i)=P(i-l) -D T (i) (D(i) -P(i-l) -D T (i)+W(i) ) _1 

P(i)=P(i-l)-K(i) -D(i) -P(i-l) (4) 
e(i|i-l)=r(i)-D(i) -S(i-l) 
S(i)=S(i-l)+K(i) -e(ili-l) 

30 In this context, K(i) is the Kalman-gain, P is the 

prediction state covariance matrix, D is the data matrix, 
which results from (2) and/or (3), W is the noise- 
covariance matrix and S is the vector of the estimated 
scatterer coefficients, which results from the 
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arrangement of the scatterers in a linear vector from the 
matrix S. r(i) is the received, sampled signal value 
(time or frequency domain) , i is the index in the time or 
frequency direction . 

5 

The methods of recursive estimation are per se known and 
have been described, for example, in S. Haykin, "Adaptive 
Filter Theory", 1 st Edition, Englewood Cliffs, New Jersey, 
Prentice Hall 1986. 

10 

It should also be noted that the RLS algorithm described 
is only mentioned as one example of a large number of 
different embodiments . 

15 After the initial estimation of the channel using 

training symbols, a maximum likelihood (ML) approach is 
selected, in which minimisation is carried out in the 
equations (2) and (3) for unknown data sequences over all 
possible data sequences and all possible arrangements of 

20 scatterers. 

A tree-search procedure can advantageously be used in 
conjunction with the channel estimation. In this context, 
starting from the channel estimated with reference to the 

25 training sequence, a pathway in a tree is built up by the 
receiver for each of the potential data sequences. A 
channel estimation is carried out with the estimation of 
the scatterers for each of these pathways, and a metric 
is calculated according to equation (2) and/or (3). The 

30 data sequence with the best metric is presented as the 
data sequence which has probably been received. Because 
of the ML approach, the metric is known as a ML-metric. 
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Instead of the metrics according to (2) and/or (3), which 
are determined in one block over the entire observation 
interval N, an incremental metric may also be used. This 
takes equation (4) into consideration as follows: 

5 

A(i)= A(i-l)+e(i|i-l) ■ (r (i) -D (i) H S (i) ) (5) 

This tree-search procedure is illustrated schematically 
in F igured for binary symbols, A,(x,...y) denotes the 

10 metric for the assumed symbols x..y, S denotes the matrix 
for the scatterers determined for the relevant pathway. 
The number of indices indicates the depth of the tree; in 
the example, up to a maximum of three. The additionally 
marked pathway characterises the best pathway selected 

15 via the metric at the moment. 

The algorithm described is a soft output algorithm, 
which, alongside the demodulated data, can also present a 
quality measure for the demodulation in the form of a 

20 metric. Accordingly, it is possible to present not only 
the data sequence determined as the most probable, but 
also less probable sequences. Processing stages, such as 
decoders connected downstream in the receiver, can 
contain additional information, which has a positive 

25 influence on the quality of reception. 

In this manner, several data sequences can continue to be 
processed in the subsequent processing stages, a decision 
about the actual sequence received being made only 
30 afterwards. 



Moreover, the method can advantageously be combined with 
a convolutional code or a block code as a single code or 
internal code of a concatenated code structure. 
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Presentation of convolutional codes and block codes in 
the form of tree structures is already known. A code acts 
on the above-mentioned tree structure in such a manner 
that not all pathways, which would be possible if the 
5 code were not taken into consideration, actually exist. 
Accordingly, when code information is included, a tree of 
this kind will not contain all pathways. 

This combination provides combined channel estimation, 
10 equalisation, demodulation and decoding, which is 

referred to as "sequential decoding". Although this 
method is already known, its use in conjunction with the 
determination of the scatterer coefficients is novel. 

15 A tree derived from the example of Figure 2 is 

illustrated in Figure 3. Comparison of the two trees 
shows that pathways determined by the code are non- 
existent . 

20 With multiple-value data symbols and/or long data 

sequences, very many pathways occur during the course of 
processing, for each of which the metrics and scatterer 
matrices as well as other auxiliary parameters for the 
algorithms must be calculated and stored. The number of 

25 pathways can be reduced in order to lessen the burden of 
calculation and memory requirement. In this context, the 
total number of pathways is limited to a maximum number, 
which depends on the available calculation capacity and 
the memory requirement of the receiver. In this context, 

30 the known metric-first, breadth-first or depth-first 
algorithms can be used. 

Known special methods for equalisation with a tree-search 
procedure have disadvantages in the context of channels 
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with long impulse responses, in which a large proportion 
of the energy of a data symbol is disposed at the end of 
the impulse response, so that this energy is not included 
in the estimation of the received symbol. In this 
5 context, the entire impulse response must either first be 
waited for with a corresponding, additional delay, or it 
must be taken into account through additional estimation 
methods with a modelling of these influences as noise. 
With the first variant, many additional pathways occur, 
10 which have to be included in the computation, even if 
they are rejected afterwards. If the method is used for 
general and unknown channels, the computations must 
always use the maximum channel impulse lengths, and the 
algorithm must therefore be designed for this in advance. 

15 

The method according to the invention does not avoid 
these disadvantages a priori. However, since the channel 
is modelled with reference to the scatterers, the maximum 
occurring delay, and therefore the dimension of the 

20 scatterer matrix, can be measured by determining the 
relevant scatterers. While this maximum length must 
always be taken into consideration in the context of the 
known methods, the method according to the invention 
allows the maximum delay of the channel to be approached 

25 in an adaptive manner, and the necessary delay in 

demodulation and decoding is adjusted accordingly. A long 
additional delay in the demodulation and coding becomes 
necessary only in special channels, in which significant 
scatterers occur with long delays. Since the geometry of 

30 the scatterers does not change abruptly, the dimension of 
the scatterer matrix can be increased adaptively if a 
scatterer with long delay occurs. Conversely, if a 
scatterer of this kind disappears, the dimension of the 
matrix can be adaptively reduced. 
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The decision can be represented in terms of a formula 
based on (2) : 



cf (0. .N-L-l) 



arg mm 



d(Q..N-L-\)d{N-L..N-L) 
S(m.k) 



N-\\\ 1 M K J2a *L 

i=0 V N m-0 k=-K 



2\ 



(6) 



In this context, L is the necessary delay. The minimum is 
determined for all possible data hypotheses d and all 
10 possible scatterers S. 



In addition to optimising the dimension of the scatterer 
matrix with reference to delay, the maximum Doppler shift 
occurring can also be optimised. 

15 

In the context of equalising and demodulating single 
carrier methods, the transmitted data can only cause ISI 
in the time direction, that is to say, data transmitted 
in the past influence data transmitted at a later time. 

20 

Because of the ICI occurring in the frequency domain when 
receiving multi-carrier signals, e.g. OFDM, a given 
carrier can be influenced by adjacent carriers both in 
the positive and also negative frequency direction. 

25 

It must also be taken into account that a cyclic 
continuation of the carriers occurs in the frequency 
domain. This cyclical continuation can be taken into 
consideration in the data matrix D, by defining the data 
30 symbols D(n-k) with a negative index occurring in 
equation (3) . 
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As in the context of considering long delays in the 
channel impulse response when processing in the time 
domain, this influence can be taken into account and 
compensated by including "future'' events, that is, data 
5 of higher frequencies, through a corresponding delay of 
the decisions. Here also, the scatterer matrix can be 
varied adaptively. 

An analogous decision for multiple carrier methods is 
10 achieved if (3) is used in (6) . 

The method described can also operate without 
initialisation based on training sequences. In this case, 
processing is initialised with default values, e.g., the 
15 matrix P from (4) is pre-defined as the unity matrix, and 
the scatterer vector S is initialised at zero. The 
algorithm will then generally converge more slowly. 
Furthermore, all possible starting configurations for the 
data sequence must be included. 

20 



